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Infrared and X-Ray Evidence for Circumstellar Grain Destruction 
by the Blast Wave of Supernova 1987A 
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John Danziger^, James M. De Buizer^, Robert D. Gehrz^, Robert P. Kirshner^, Richard 
McCray^, Sangwook Park^, Elisha F. Polomski^, and Charles E. Woodward^ 

ABSTRACT 

Multiwavelength observations of supernova remnant (SNR) 1987A show that 
its morphology and luminosity are rapidly changing at X-ray, optical, infrared, 
and radio wavelengths as the blast wave from the explosion expands into the 
circumstellar equatorial ring, produced by mass loss from the progenitor star. 
The observed infrared (IR) radiation arises from the interaction of dust grains 
that formed in mass outflow with the soft X-ray emitting plasma component 
of the shocked gas. Spitzer IRS spectra at 5 - 30 fim taken on day 6190 since 
the explosion show that the emission arises from ~ 1.1 x 10~^ Mq of silicate 
grains radiating at a temperature of ~ 180±^5 K. Subsequent observations on 
day 7137 show that the IR flux had increased by a factor of 2 while maintaining 
an almost identical spectral shape. The observed IR-to- X-ray flux ratio {IRX) 
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is consistent with that of a dusty plasma with standard Large Magellanic Cloud 
dust abundances. IRX has decreased by a factor of ~ 2 between days 6190 and 
7137, providing the first direct observation of the ongoing destruction of dust in 
an expanding SN blast wave on dynamic time scales. Detailed models consistent 
with the observed dust temperature, the ionization timescale of the soft X-ray 
emission component, and the evolution of IRX suggest that the radiating silicate 
grains are immersed in a 3.5 x 10^ K plasma with a density of (0.3 — 1) x 10^ cm~^, 
and have a size distribution that is confined to a narrow range of radii between 
0.023 and 0.22 /im. Smaller grains may have been evaporated by the initial UV 
flash from the supernova. 

Subject headings: ISM: supernova remnants - ISM: individual (SNR 1987A) - 
ISM: interstellar dust - Infrared: general - X-rays: general 



INTRODUCTION 



On February 23, 1987, Supernova 1987A (SN 1987A), the brightest supernova since 
Kepler's SN in 1604, exploded in the Large Magellanic Cloud (LMC). About ten years 
thereafter, the energy output from the supernova became dominated by the interaction of 
its blast wave with the inner equatorial ring (ER), a dense ring of gas located at a distance 
of about 0.7 lyr from the center of the explosion, believed to be produced by mass loss 
from the progenitor star. The ER is being repeatedly observed at optical wavelengths with 
the Hubble Space Telescope (HST) (IPun et al.ll2002l ). at X-ray energies with the Chandra 



X-ray Observatory (jPark et al.ll2006l. 120071). at radio fr equencies with the Australian Tele 
scope Compact Arr ay (ATCA) (Manchester et allboosl ). and in the mid-IR w ith the Gemini 
Sout h observatorv JSouchet et all booi l2006h and the Spitzer observatory fiBouchet et al.l 



2004 : iDwek fc Arendtll2007l ). The morphological changes in its appearance in these different 
wavelength regimes (presented in Figure 6 of McCray 2007) reveal the progressive interaction 
of the SN blast wave with the ER. The interaction regions appear as hot spots in the HST 
images, representing the shocked regions of finger-like protrusions that were generated by 
Rayleigh- Taylor instabilities in the interaction of the wind from the progenitor star with the 
ER (see Figure 2 in McCray 2007). 

The X-ray emission is thermal emission from the very hot plasma, and consists of two 
main characteristic components: a hard {kT ^ 2 keV) component representing a fast shock 
propagating into a low density medium, and a soft {kT ^ 0.3 keV) component representing 
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a dec elerated shock propagating into the denser protrusions flPark et al.l 120071 : iZhekov et al. 
20061 ). The optical emission arises from the gas that is shocked by the blast wave trans- 
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mitted through the dense protrusions in the ER (jPun et aLll2002l ). and t he radio emission 



is syn chrotron radiation from electrons accelerated by the reverse shock (IManchester et al. 



20051 ). The mid-IR emission spectrum is comprised of line and continuum emission. The 
lines most probably originate from the optically bright dense knots. The continuum that 
dominates the spectrum is thermal emission from dust that was formed in the post main 
sequence wind of the progenitor star before it exploded. This dust could either be located in 
the shocked X-ray emitting plasma and heated by electronic collisions or in the optical knots 



and radiatively heat ed by the shocks giving rise to the optical emission (jPolomski et al.ll2004 



Bouchet et al.ll2006l ). Dust has also formed in the ejecta of the supernova a bout 530 days af- 



ter the explosion, as evidenced by optical and IR observation s of SN1987A (ILucy et al.lll991 



Gehrz fc NevlllQQol : iMoselev et al1ll989l : Iwooden et aljll993h . but its current contribution to 
the total mid-IR emission is negligible (IBouchet et al.ll2004l ). 



Because the Gemini 11.7 /im image correlates well with both the X-ray and optical 
emission, th e possibility that the dust is radiatively heated in the knots was considered 
in detail by iBouchet et al.l (120061 ). Estimated dust temperatures of ~ 125 K fell short of 
the observed value of ~ 180 K but given the uncertainties in the model parameters, this 
discrepancy could not firmly rule out this possible scenario for the location and heating 
mechanism of the dust. However, the combined IR and X-ray observations can be more 
readily explained if the dust resides in the shock ed regions of the fi nger- like protrusions 
that give rise to the observed soft X-ray emission (jZhekov et al.ll2006l ). We will adopt this 
scenario as our working hypothesis, and derive a self-consistent model for the composition, 
abundance, and size distribution of the dust to explain the evolutionary changes in the IR 
and X-ray fiuxes resulting from the collisional heating and the destruction of the dust grains 
by the ambient plasma. 

We first review the basic physical principles that determine the temperature of coUi- 
sionally heated dust, and describe how the IR emission can be used to probe the physical 
condition of the X-ray emitting plasma. We also discuss what information can be derived 
from the comparison of the IR and X-ray fiuxes from the gas (§2.1). The physics of dust 
particles in a hot gas is discussed in more detail by Dwek (1987) and Dwek & Arendt (1992). 
In §3 we present a simple analytical model for the evolution of the grain size distribution 
and total dust mass in the gas that is swept up by an expanding SN blast wave. In §4 
we present Spitzer low resolution 5 — 30 /im IRS spectra obtained on days 6190 and 7137 
after the explosion. The IR spectra are used to derive the temperature and composition 
of the shock heated dust. We use IR and X-ray observations of the SN to constrain the 
grain size distribution and the time at which the SN blast wave first crashed into its dusty 
surroundings. The results of our paper are summarized in §5. 
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2. The Infrared Diagnostics of a Dusty X-ray Plasma 



The morphological similarity between the X-ray and mid-IR images of SN 1987A sug- 
gests that the IR emission arises from dust that is coUisionally heated by the X-ray emitting 
gas. Simple arguments presented below show that, under certain conditions, the IR luminos- 
ity and spectrum of a dusty plasma can be used as a diagnostic for the physical conditions of 
th e gas and the det ails of the gas-gra i n inte ractions. Details of the arguments can be found 



m 



Dwekl (Il987f ) and lPwek fc Arendtl fll992[ ). 



2.1. The Dust Temperature as a Diagnostic of Electron Density 

The collisional heating rate, 7i(erg s~^), of a dust grain of radius a embedded in a hot 
plasma is given by: 

n = Txa^ ^rij Vj £,j (1) 
i 

where Uj is the number density of the j-th plasma constituent, Vj, its thermal velocity, and 
£j its thermally- averaged energy deposition in the dust. In all the following, we will assume 
that the ion and electron temperatures are equal. Then Vf. ^ Vion, and the dust heating rate 
is dominated by electronic collisions. 

Let Edep be the thermally-averaged energy deposited by electrons in the solid. If most 
electrons are stopped in the dust then Edep is, on average, equal to the thermal energy of 
the electrons, that is, Edep oc Te, where Tg is the electron temperature. On the other hand, 
if most incident electrons go entirely through the grains, then Edep is proportional to the 
electron stopping power in the solid, defined as p~^{dE/dx). At the energies of i nterest here 



the e lectronic stopping power has an energy dependence of (dE/dx) ~ E ( Iskef et al. 



1983|), or dE ~ E~^/^ dx so that Edep ~ Te 
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a. 



The functional dependence of the grain heating rate on gas density and temperature is 
then given by: 

H. ~ Tie T^^"^ when electrons are stopped in the grain (2) 

~ 0? He when electrons go through the grain 

1/2 

where we used the fact that Ve ^ Te . 

The radiative cooling rate, £(erg s~^), of the dust grain with temperature Td by IR 
emission is given by: 

C = na' an iQ) (3) 
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3 A7i4+/3 

where a is the Stefan-Boltzmann constant, and (Q) oc aT^ is the Planck-averaged value of 
the dust emissivity, Q(A) oc A"'^, where the value of emissivity index, is ~ 1 — 2. 

In equilibrium, H = C, and the dust temperature dependence on plasma density and 
temperature can be written as: 

Td ~ \ ~) '^t'^^'^ when electrons are stopped in the grain (4) 

~ nl when electrons go through the grain 

where 7 = 1/(4 + These simple arguments show that when the gas temperature is 
sufficiently high, and the grain size sufficiently small so that most electrons go through the 
grain, the dust temperature depends only on the plasma density. 




10* 10^ lO** 10^ 10" 10* 10^ 10*^ 10^ 10' 

Te(K) T,(K) 



Fig. 1. — Contour plot of the equilibrium temperature of 0.01 /xm (left) and 0.10 /.im (right) sihcate grains 
as a function of electron density and temperature. Above temperatures of ~ 5 x 10^ K (~ 3 x 10'' K) the 
0.01 /im (0.1 /im) grains become transparent to the incident electrons, and the dust temperature is only a 
function of electron density. At lower temperatures, different combination of plasma density and temperature 
can heat the dust to the same temperature. We note here that a similar figure (Figure 15 in Bouchct ct 
al. 2006) was mislabeled, and actually corresponds to contour levels of silicate dust temperature for grain 
radius of a = 0.0030 /im. 

Figure [H depicts contour levels of the dust temperature as a function of electron density 
and temperature for 0.01 //m and 0.10 [ira. silicate grains. The figure shows that above a 
certain gas temperature, its value depending on the grain radius, most of the electrons go 
through the grain and the dust temperature is essentially determined by the electron density. 
Under these conditions, the IR spectrum from the collisionally-heated grains becomes an 
excellent diagnostic of the density of the X-ray emitting gas. However, when the grains 
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are large enough to stop the electrons, the dust temperature only constrains the possible 
combination of plasma temperature and density. 



2.2. The Stochastic Heating of Grains by Electronic Collisions 



When dust grains are sufficiently small, a single electronic collision can deposit an 
amount of energy in the dust that is significantly larger that its enthalpy, causing a surge in 
dust temperature. If additionally, the time interval between successive electronic collisions 
i s larger than the dust cooling tim e, the grain temperature will be fiuctuating with time 
(iDwekl 119871 : iDwek fc Arendtlll992l ). Figure [2] depicts a simulation of the stochastic heating 
of 0.003, 0.01, and 0.03 /im silicate grains immersed in a hot X-ray emitting gas characterized 
a temperature Tg = 10^ K, and an electron density Ue = 1 cm~^. The left column shows 
the temperature fluctuations as a function of time, and the right column the histogram of 
the grain temperature. As the grain size increases, the fluctuations get smaller, and the 
histogram becomes strongly peaked around the equihbrium dust temperature of ~ 38 K, in 
this example. 



2.3. The Infrared to X-ray Flux Ratio 

Another important diagnostic of a dusty plasma is IRX, deflned as the ratio of the IR 
to X-ray fluxes emitted by the gas (Dwek et al. 1987). If the dust is collisionally-heated 
by the gas then the total IR flux, F^, emitted from a gas volume V is proportional to 
He rid ^d{Tg) V, where rid is the number density of dust particles, and Ad{Tg) is the cooling 
function (units of erg cm^ s~^) of the gas via gas-grain collisions. The total X-ray flux. Fx, 
from the same volume is proportional to Ag{Tg) V, where Ag{Tg) is the cooling function 
of the gas via atomic processes. Thus 

lUX^C^)^ (5) 

Both cooling functions represent the energy losses through coUisional processes, characterized 
by {av E) summed over all interactions in the plasma, where a is the cross section, v is the 
relative velocity of colliding species, and E is the energy lost in the process. 

For a given dust-to-gas mass ratio, that is, a flxed (n^/rie) ratio, IRX depends only on 
plasma temperature. Figure [3] (left panel) shows the behavior of the atomic cooling function 
of a gas of solar composition as a function of gas temperature. Also shown in the flgure 
is the gas cooling function via gas-grain collisions for a gas with a dust-to-gas mass ratio 
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Fig. 2. — The stochastic heating of silicate grains in a hot X-ray emitting gas characterized by a temperature 
of Tg = 10*' K, and electron density rie = 1 cm^'^ for dust grains of different radii. Left column: The 
temperature fluctuations as a function of time. Right column: The histogram of the fluctuations. As 
the grain size increases, the fluctuations get smaller, and the probability distribution of dust temperatures 
becomes strongly peaked around the equilibrium temperature of 38 K. 
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Fig. 3. — Left panel: The cooling function of a dusty plasma via atomic processes (thick solid line) and 
via gas grain collisions. Calculations were performed assuming a single-sized population of 0.02, 0.1 and 
0.20 /im grains (solid line) with a dust-to-gas mass ratio of 0.0062, which is the value in the local ISM of 




the bare silicate-graphite+PAH dust model of IZubko et al.l ([200J). The gas cooling rate per unit volume for 
both processes is given by L = A{T). Right panel: The value of IRX for single-sized dust populations 
with radii of 0.02, 0.1, and 0.2 jum with the same dust-to-gas mass ratio as the figure on the left. 



Zd = 0.0062 fIZubko et al.ll2004l ). and single-sized dust populations with radii of 0.02, 0.1, 
and 0.2 /xm. The right panel of the figure presents the value of IRX for the same conditions. 
The figure shows that for soft X-rays {kTe ~ 0.3 keV, Tg ~ 3.5 x 10^ K) this ratio varies 
between ~ 3 and 20, depending on grain size. Each plasma temperature will have a different 
range of values, depending on the grain size distribution. Any deviation from these values 
will suggest that Zd is either depleted or overabundant with respect to the reference value 
adopted in the calculations. 



2.4. The Plasma Ionization Timescale 

The cooling rate of a plasma depends on the ionization state of its constituent ions 
which may not have evolved to equilibrium conditions. The ionization state of the gas is 
characterized by the ionization timescale, jF(units of cm~^ s), defined as = rig x te, where 
te is the age of the shocked gas. In a fully ionized plasma with ~ Uion, T also measures 
the fluence of ions incident on the dust. When sputtering is the dominant grain destruction 
mechanism, T will directly determine the total mass of dust that is returned to the gas. 

In summary, the key parameters: plasma temperature, density, ionization timescale, and 
X-ray fluxes, and the dust temperature, grain size distribution, composition, and IR fluxes 
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are closely interrelated so that knowledge of some parameters will constrain the others. 



3. The Evolution of the Grain Size Distribution and Dust Mass 

Consider the propagation of a shock into a dusty medium with a constant number density 
and a fixed dust-to-gas mass ratio Z^. The shocked gas can be regarded as a reservoir that is 
continuously being filled with gas and pre-existent circumstellar dust by the expanding blast 
wave. If dust grains were not destroyed, the postshock gas would maintain a constant value 
as the mass of shocked dust and gas evolve proportionally in time, with a functional 
dependence that depends on the geometry of the medium into which the blast wave is 
expanding. In the case where dust particles are destroyed by sputtering, the grain size 
distribution, and the dust-to-gas mass ratio in the shocked gas will evolve with time. 



3.1. General Equations 

A grain of radius that is swept up by the shock at some time to will at time t be 
eroded to a radius a{t) given by 

a = a{t) = «o + ^ (J^^ dt' (6) 

where 

da V — > 

-~ J] n,t;,F, <0 (7) 
j 

where Yj is the thermally-averaged sputtering yield of the dust by the j-th gas constituent. 

Dust particles swept up by a high velocity shock will move ballistically through the 
shock front and acquire a velocity relative to the shocked gas. The sputtering yield needs 
then to be averaged over a Maxwellian distribution of ve locities that is displa ced by the 



relative gas-grain motion from its origin in velocity space (IDwek fc Arendtlll992l ). Figure H] 
shows the temperature dependence of the sputtering rate of silicate dust grains, calculated 
using sputtering yield parameters given by Nozawa et al. (2006), moving with velocity Vgr 
= 0, 50, 100, 200, and 500 km s~^ through a hot gas of solar composition. In contrast to the 
heating of grains, their erosion by thermal and kinetic sputtering is entirely done by the ionic 
constituents of the gas. For dust grains with velocities > 500 km s~^ and gas temperatures 
above ~ 10^ K the sputtering rate is approximately constant and given by: 

^ ^ -5 X 10"*^?7,H(cm"^) /im yr"^ (8) 
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Fig. 4. — The absolute value of the sputtering rate (normalized to the H-number density) of silicate 
(Mg2Si04) dust grains moving through a hot gas of solar composition as a function of gas temperature. The 
curves are marked by the grain velocity (in km s^^). 



Equation ([H]) ignores the possible destruction of grains by evaporative grain-grain collisions 
that occurs in slower (< 100 km s~^) shocks by the acceleration of charged grain with MHD 
turbulence (lYan et al.ll2004l ). 

If the shocked gas maintains a constant composition and density then a dust grain of 
initial radius that is swept up by the shock at some time t' will at time t have a radius a 
given by: 

(9) 



Equation ([9]) can be written in dimensionless form as: 

(t - t') 



where ^ = a/a^ax, Co = ^o/amax, V = t/rm.ax, v' = i'/Tmax, and 

'Tmax — Oimax \dci/dt\ 



(10) 



(11) 



is the sputtering lifetime of the largest grain in the injected size distribution, which (using 
eq. [HI) is numerically given by: 



2 X 10' 



(12) 
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7.3 X 10 
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nHicm ^1 



Let Ud be the total number density of dust grains in the preshocked gas, and nd{ao) daQ, 
the number density of grains with radii between cq and + daQ. We assume that the grains 
have a size distribution in the preshocked gas given by: 

nd= nd{ao)dao = nd / f{ao)dao (13) 
Jo Jo 

where /(cio) is the normahzed size distribution. If the grain size distribution extends over a 

hmited range of radii, amin < cio < amax, then /(ao) = for any ao < a™„ or ao > amax- 

Dust grains are continuously injected into the shocked gas by the expanding SN blast 
wave. The total number of shocked grains with radii a in the {a, a + da\ interval at time t, 
Nd{a,t)da, is equal to the number of all dust particles of initial radius Qq that were swept 
up at time t' (0 ^ t' ^ t) and sputtered during the time interval t' — t to radius a given by 
eq. (Q. If V(t) is the growth rate of the volume of the shocked gas, then Nd{a,t) can be 
written as: ^ 

Nd{a,t) = nd I V{t')f{ao)dt' (14) 



The lower limit of the integral, t = 0, corresponds to the time when the blast wave first 
encounters the dusty medium. 

The total mass of shocked dust at any given time t is given by: 



Md{t) = I md{a)Ndia,t)da (15) 

„ — \da/dt\t 



where md{a) = Anpa /?> is the mass of a dust grain of radius a. 



Equation can be written in dimensionless form: 



Nd{i, V) = rid 



rv 


dV{r]'Tmax) 


Jo 


drj' 



f[{^ + n-r^')a,,ax]dr,' (16) 



This integral is a convolution of the form: g{ri') * f{rio — f]'), which can be numerically 
evaluated for arbitrary functions using Fourier transforms. 



3.2. A Simple Analytical Solution 



An analytical solution can be derived for a pre-shocked grain size distribution given by 
a power law in grain radius, and a power law time dependence of V. We write the grain size 
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distribution as: 



/(«o) 



otherwise 



(17) 



fc+i 



where C = {k — l)/(a. 
The time dependence of V{t') can be written as: 



V{t') = Vo 



is the normalization constant. 



t' 



r„ 



where Vq is a proportionality constant, and a = 2 for a spherical blast wave expanding 
into a uniform interstellar medium (ISM), and a = if the blast wave expands into a 
one-dimensional "finger-like" protrusion. 

The total number density of grains in the {a, a + da} radius interval is then given by: 



Nd{a,t) = UdVo 



/(ao) dt' 



(19) 



'0 \ rnax . 

Using eq. ([9]) to change variables from t' to ao, eq. ( fT9l) can be rewritten as: 



Nd{a,t) = Nd 



da 


~1 rdup 






[( 


dt 


(^low 





t a 
+ 



max ""niax 



ao 



(3g ''^ dao 



(20) 



where A^^^ = Vq. 



The time dependence of Nd{a, t) is contained in the limits on the integral over the grain 
size distribution. If the radius a is within the range of the injected grain size distribution. 



that is, 



< a < ar, 



then aiow = cl, since only grains with radii larger than a could have 



contributed to Nd{a,t). If the radius a is smaller than amin, then Nd{a,t) is non-zero only 
if a -|- \da/dt\t exceeds amin, and aiow = CLmin- In other words, the most recent injection of 
grains that could have contributed to Nd{a, t) occurred at time t — At, where At is the time 
required to reduce the grain radius from amin to a. The largest grains that could have been 
sputtered to radius a during the time t is equal to a + \da/dt\t. However, the largest grain 
size cannot exceed amax, so the upper limit on the integral, aup, is determined by the smaller 
of these two quantities. To summarize: 



max{amm) o,} 



^up mm amax i O' 



(21) 



da 



dt 
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For a spherical shock wave expanding into a one-dimensional protrusion [a = 0) the 
solution to eq. fl20|) is given by: 



Nd{a,t)^=o = Nd 



da 



dt 



C 



{k 



^\ ["'low 



-k+1 



up 



(22) 



At early times, when t < a\da/dt\ ^ ^ Tmax, dup = cl, and aiow ~ cl, and the solution to eq. 
(!22|) becomes: 

Ndia,t)^=o = NdtCa-'' . (23) 

At late times, when t > Tmax, o-up = dmax, and the asymptotic solution of ( l22i) is: 

-1 



Nd{a,t> 



da 



dt 
da 



C 



dt 



{k- 
C 



^max 



min max 



] a > amin (24) 

] = constant a < amin 
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Fig. 5. — Evolution of the grain size distribution with time, measured in units of Tmax, the sputtering 
hfetime of the largest grain with radius amax hi the size distribution. Calculations were performed for a 
grain size distribution characterized by an a~^'^ power law in grain radii between 10 and 100 A. The 
grain destruction rate, da/dt, was taken to be 0.14 A d~^, for an assumed density of 10^ cm^^. Bold lines 
are labeled by t/Tmax- Left: A spherical blast wave expanding into a one-dimensional protrusion (a = 0). 
Right: A spherical blast wave expanding into a uniform ISM (a = 2). 



For a shock wave expanding into a homogeneous medium (a=2) the solution is given 



by: 



Nd{a,t) 



a=2 



Nd 



da 






dt 



{k 



- a~''+^ 

ow up 



) (25) 
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{k-2) \t, 



t a 
+ 



-fc+2 _ „-fc+2' 



+ 



' max ^max 



-/c+3 _ „-/c+3' 
"'low "'up 



{k-3) 



At late times, when t ^ Tmax, the first term dominates, and the asymptotic solution of eq. 
f l25|) increases with time as t^: 



Nd{a,t > Tmax)a=2 = Nd 



da 



dt 



C 



(k-l) Vr, 



t 



where 



= [(^min^ - ^mat^] = COnstaut 



(I ^ Q'min 
Qj O'TD-in. 



(26) 



(27) 



Figure \5\ depicts the grain size distribution for different epochs. Select epochs, labeled 
by the dimensionless quantity t/Tmax, are represented by bold lines. Calculations were per- 
formed for an initial grain size distribution characterized by an ~ a~^'^ power law in grain 
radius between 10 and 100 A. The grain destruction rate, \da/dt\, was taken to be 0.14 A 
d~^, calculated for the sputtering rate of silicate grains in a hot gas with a temperature and 
density of ~ 10^~^ K, and 1000 cm~^, respectively 

The figure illustrates the dependence of the evolution of the grain size distribution on 
the geometry of the ISM into which the blast wave is expanding. For a one- dimensional 
protrusion (a = 0), the figure (left) shows a clear convergence of the size distribution to a 
fixed functional form and total number of grains for t/Tmax ^ 1- As the shock wave expands, 
the thickness of the shell of swept up dust increases with time. However, because of the 
finite grain lifetime, its thickness cannot exceed a value of ARgh ~ Vgh Tmax, where Vsh is the 
shock velocity. So the grain size distribution and total mass reaches a steady state limit. 
When the blast wave expands into a uniform medium, the shell of shocked dust reaches the 
same steady state thickness ARgh- However, since the surface of the shell increases as R'^f^, 
where Rgh is the radius of the blast wave, the mass of shocked gas will continue to increase. 
This is clearly depicted in the right panel of the figure, which shows that Nd{a,t) reaches a 
steady state, but continues to increase with time as t^. 

If grains were not sputtered in the shocked gas, then N^{a,t) da, the total number of 
dust grains in the {a, a + da} radius interval that are swept up by the shock at time t would 
be: 

c ^ t \°''^^ 

N^{a, t) da = Nd Tmax T da (28) 

(a + 1) \TmaxJ 
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Fig. 6. — Evolution of the dust mass, M^, the dust mass if grains were not sputtered, M^, and the fraction 
of the surviving dust, M^/M^ with time, measured in units of Tmax, the sputtering hfetime of the largest 
grain in the size distribution [Trnax{d) = 7.3 x 10'' a,„Q2:(Mm)/n/f (cm^'^)]. Dust and gas parameters are 
identical to those used in Figure [5l The figure shows that the fractional change in Mj,/M^ between two 
epochs constrains the grain size distribution and the density of the X-ray emitting plasma. The fraction 
Md/M^ is proportional to IRX, the IR-to-X-ray flux ratio of the shocked gas. Left: A spherical blast 
wave expanding into a onc-dimcnsional protrusion (a = 0). Right: A spherical blast wave expanding into 
a uniform ISM {a = 2). 

Figure M shows the evolution of dust mass for a spherical blast wave expanding into a 
one-dimensional protrusion [a = 0; left), and into a uniform ISM [a = 2; right) as a function 
of t/Tmax- As explained above, when a = 0, the mass of shocked dust, Md, reaches a constant 
limit for t ^ Tmax, whereas for a = 2 the mass of the shocked dust will increase as t"^. If 
grains were not destroyed, the mass of swept up dust, M° would increase as t for a = 0, 
and as for a = 2. The figure also shows the evolution of the mass fraction of surviving 
dust grains, M^/M^. This mass fraction is proportional to the dust-to-gas mass ratio in the 
shocked gas, and for a constant gas temperature and density, it is also proportional to IRX, 
the IR-to-X-ray flux ratio in the shocked gas. The figure shows that the fractional change 
in M^/M^ between two epochs constrains the value of Tmax given in eq. (fT2l) which in turn 
depends on the grain size distribution and the density of the X-ray emitting plasma. As 
a reminder, Tmax{d) = 7000 amax{A)/nH{cm~^). For example, given a plasma density, the 
value of Tmax will depend only on amax, the maximum grain radius. A small value of amax 
will imply a low value for Tmax, so that large changes in Md/M^ occur over very short time 
scales. Conversely, large values of amax and Tmax, will cause changes in Md/M^ to occur over 
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very long time scales. 



4. Spitzer Infrared Observations of SNR 1987A 
4.1. The Evolution of the IR spectrum 



Figure [7] shows the 5 — 30 /im low resolution spectra of SNR 1987A taken on February 
4, 2004 (day 6190 since the explosion), and or i September 8, 200 6 (day 7137 since the 
explosion) with t he Infrared Spectrograph (IRS) ( Houck et al. 2004a| jbl) on board the Spitzer 
Space Telescope (jWerner et al.ll2004l : iGehrz et al.ll2007l ). Analysis of the spectrum taken on 
day 6190 revealed that the IR emission originated from ~ 1.1 x 10~^ Mq of silicate grains 
radiating at a temperature of ~ ISOl^g K (Bouchet et al. 2006). These circumstellar grains 
were formed in the quiescent outflow of the progenitor star before it exploded. The total IR 
flux on day 6190 was 5.1 x 10^^^ erg cm~^ s~^ (Bouchet et al. 2006), and increased after 
947 days (day 7137) to 10.0 x 10"^^ erg cm""^ s^^. The right panel of flgure [7] presents a 
comparison between the two spectra, both normalized to the same 10 fim intensity. The 
flgure shows that the spectra are essentially identical, implying that the dust composition 
and temperature remained unchanged during the two observing periods. The lower curve in 
the flgure is the ratio between the two spectra, emphasizing their similarity. The IR intensity 
increased by a factor of 2 between the two epochs. 



4.2. Determining the Grain Size Distribution 

The size distribution of coUisionally-heated dust grains is constrained by the combi- 
nations of gas temperature and density that can give rise to the range of observed dust 
temperatures. The temperature of the gas giving rise to the soft X-ra y componen t can be 



derived from models, and is equal to ~ 0.3 keV (Te = 3.5 x 10^ K; iPark et all (120051 ) ). 
narrowing down the range of viable plasma densities and grain sizes. 

Figure [8] depicts contours of the dust temperature as function of gas density and grain 
size for the given electron temperature, Tg = 3.5 x 10^ K. The range of observed dust 
temperature falls between 165 and 200 K, and the flgure shows the different combination of 
grain size and gas density that can give rise to this narrow range of dust temperature. All 
grain sizes are viable, provided that the gas density has the right value to heat the dust to 
observed range of temperatures. However, the range of viable grain sizes can be narrowed 
down by using the constraints on the ionization timescale of the plasma. 
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Fig. 7. — Left: The evolution of the IR spectrum of SN1987A from April 2, 2004 (day 6190 since the 
explosion) to September 8, 2006 (day 7137) taken with the Spitzer IRS (Bouchet et al. 2007; Arendt et 
al. 2007). Right: The smoothed spectra for days 6190 (solid line) and day 7137 (dashed line), normalized 
to the same brightness. The lower curve shows the ratio between the two spectra, with the horizontal line 
being the mean value. The figure shows that the dust spectrum increased by a factor of two between the two 
epoch, retaining essentially an identical spectrum corresponding to silicate grains radiating at an equilibrium 
temperature of 1801^" K. 



The ionization timescale derived from modeling the soft X-ray spectr a taken on days 



6914, 7095, and 7271 is given by JF = ngte ^ 10'' cm ^ d fjPark et aLl 120071 ). The ionization 



time, te, is constrained by the time to when the SN blast wave first encountered the ER. 
Observationally, we can associate to with the appearance of t he first hot spot in the HST 
image from April 1997, about 3700 days after the explosion (jPun et al.l l2002l ). The soft 
X-ray light curve shows that the rise could have occurred between days 3700 and 6000. The 
first epoch corresponds to the first appearance of the optical knots, and the latter epoch 
corresponds to the time when the flux from the soft X-ray componen t {kT ~ 0.3 keV) 
exceeded tha t from the hard cora ponent {kT ~ 2 keV) (jPark et al.ll2005l ). From the mid-IR 
light curves (IBouchet et al.ll2006l ). the energy output from the SN became ER dominated 
around day 4000. Adopting days 4000 to 6000 as a reasonable estimate for to gives a range 
of possible ionization times of te ~ 7000 — to ~ 1000 — 3000 d. Using the constraints on JF, 
the corresponding limits on the electron densities are: Ue ~ 10^ — 3 x 10^ cm~^. 

These densities are high enough that even small dust grains with radii ~ 10 A will be 
collisionally heated to their equilibrium dust temperature. Furthermore, for a gas temper- 
ature of Te !^ 3.5 X 10^ K, an equilibrium temperature of ~ 180 K can only be reached at 
these high densities if the soft X-ray electrons are stopped in the grains (see Figure Hj). The 
narrow range of grain temperatures then suggests that the grain size distribution should 
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Fig. 8. — Contours of the dust temperature as function of gas density and grain size for the given electron 
temperature, Tg, of the soft X-ray component. For the range of electron densities constrained by the plasma 
ionization timescale (indicated by the dashed horizontal lines), the observed range of dust temperature 
(indicated by the thick contours) limits the grain size distribution to be between ~ 0.023 — 0.22 /im. 



have a narrow range as well, since ~ a~'^ (see eq. HI). Figure E] shows that the constraint 
on the electron density, ~ (0.3 — 1) x 10^ cm~^, limits the range of viable grain sizes that 
can be heated up to ^ 165 — 200 K to be between 0.02 3 and 0.22 /im. This range is 
narrower than that adopted by IWeingartner fc Draind (120011 ) to model the size distribution 
of LMC silicate dust (see §4.7). 

The smaller upper limit on the size distribution of the silicates in the ER may be limited 
by the nucleation time scale in the outflow of the SN1987A progenitor. The higher lower 
limit on the silicate size distribution in the ER may be the result of evaporation by the initial 
UV flash (see §4.3). E.ST images of the ER show that it is located at a distance of ~ 0.7 lyr 
(6.6 X 10^^ cm) from the SN. At this distance sr aall dust particles ca n be evaporated by 
the initial UV flash that emanated from the SN (IFischera et al.ll2002l ). Their calculations 
suggest that silicate dust particles with radius less than ~ 0.02 /im will be evaporated by 
the flash. A population of silicate grains with a a""^'^ power law distribution in grain radius 
extending from 10 A to 0.2 /im will loose about 30% of its mass. 
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4.3. The Value of IRX: Constraining the Dust Abundance in the ER 



A comparison of the IR and X-ray fluxes provides strong constraints on the dust abun- 
dance in the sho cked gas. X-ray fl uxes taken between days 6157 and 7271 with the Chandra 
X-ray telescope (jPark et all 120071 ) were interpolated for days 6190 and 7137 of the Spitzer 
observations. The total X-ray flux on day 6190, corrected for an extinction H-column density 
of Nh = 2.3 5 X 10^^ crn"^, is 2.1 x 10~^^ erg cm~^ s"^, half of it radiated by the slow shock 
component (jPark et al.ll2005l ). The IR emission originates from the slow shock component 
which is penetrating the denser regions of the ER. This component comprises half of the 
observed X-ray flux. The resulting value of IRX on day 6190 is therefore 4.9 ± 1.1. 

The th eoretical value for IRX in a gas with LMC ISM abundances, taken here to be 0.6 
times solar (jWelty et al.lll999l ). ranges from about 2 to 12 for soft X-rays with Te ~ 3 x 10^ K. 
The dust abundance in the ER is therefore consistent with LMC abundances. Since the 
silicon abundance in the ER should not have been altered by stellar nucleosynthesis, this 
agreement suggests eflicient condensation of silicate grains in the presupernova outflow. The 
dust abundance on day 6190 is a lower limit on the original pre-SN value, since some of the 
dust may have been evaporated by the initial UV flash from the SN. 



4.4. The Evolution of IRX 

4.4-1- Evidence for Ongoing Grain Destruction by the SN Blast Wave 

In §3 we consider the propagation of a shock into a dusty medium with a constant 
number density and a flxed dust-to-gas mass ratio. If dust grains were not destroyed, then 
postshock gas would maintain a constant dust-to-gas mass ratio, and the X-ray and IR fluxes 
from the shocked gas will evolve proportionally in time, that is, the value of /i?Awill remain 
constant. Any evolution in the value of IRX should therefore suggest a breakdown in the 
assumptions of the model. 

Observational evidence, summarized in Table 1, show that the IR flux increased by 
only a factor of ~ 2 from day 6190 to 7137. In comparison, the extinction-corrected 0.50- 
2.0 keV flux increased by a factor of ^ 3 during the same time period to a value of ~ 



6.4 X 10 erg cm ^ s ^ fjPark et al.l 120071 ). The fractional contribution of the soft X- 



ray component increased from 0.5 to 0.6, with no signiflcant change in gas temperature 
{kT ~ 0.3 keV). All the increase in the soft X-ray flux can therefore be ascribed to an increase 
in the volume of the dense (ng ~ lO'^ cm~^) component of the ER that was shocked by the 
SN blast wave. The evolution in the X-ray and IR fluxes and the resulting value of IRX are 
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summ arized in Table 1. A similar evolutionary trend was reported by lBouchet et al.l (120041 . 



20061 ) although absolute values of IRX differ from those reported here because of differences 



in the X-ray energy bandpasses and dust models used in the calculations. 





Fig. 9. — Evolution of the ratio IRX{t2)/IRX{ti) as a function of the time difference ti — to, The time 
to is the time since the explosion when the SN blast wave first encountered the dusty equatorial ring (ER), 
taken here to be the independent variable. The times ti and t2 correspond, respectively, to days 6190 and 
7137 since the explosion, the two epochs of near-simultaneous Spitzer and Chandra observations of the ER. 
The dashed horizontal line depicts the nominal value of the I RX {t2) / 1 RX [ti) ratio which is 0.53 ±0.16 (see 
Table 1). The curves are labeled by amax, the maximum grain size (in A) of the distribution. The minimum 
grain size was taken to be 10 A in all cases. The red curves correspond to a narrow grain size distribution of 
the ER with {a,„i„, amax} = {0.023, 0.22} /im. Figures are also labeled by the value of Ugas, the density of 
the shocked gas. Left column: A spherical blast wave expanding into a one-dimensional protrusion (a = 0). 
Right column: A spherical blast wave expanding into a uniform ISM (a = 2). 
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If grains were not destroyed, we would expect the IR intensity to increase by a similar 
factor. The smaller increase in the IR flux, that is, the decline in IRX , is a strong indicator 
that we are for the first time witnessing the actual destruction of dust in a shock on a 
dynamical timescale! If the dust composition and size distribution is uniform throughout the 
region of the ER that has been swept up by the shock, then IRX is directly proportional to 
the dust-to-gas mass ratio, Zd, or to Md/M^, the ratio between the actual mass of dust in the 
shocked gas, and the dust mass if grains were not destroyed. The magnitude of the decrease 
in IRX between the two epochs is strongly determined by the following factors: (1) the 
grain size distribution in the preshocked gas of the ER; (2) the density of the X-ray emitting 
gas and the grain destruction efficiency, which determine the rate of grain destruction; and 
(3) the total time the dust is exposed to the flux of ions, which is determined by the time 
the SN blast wave first crashed into the ER. The latter factor determines the total mass of 
dust that is returned to the gas. 

The dependence of the evolution of IRX on the grain size distribution is somewhat 
subtle. If the preshocked grain sizes are too large, then the fractional mass of the dust that 
could be destroyed during the time interval of 947 days will be too small to account for the 
observed decrease in the value of IRX. Conversely, if the grain sizes were too small, most of 
the dust mass would be destroyed, giving rise to a significantly larger than observed decrease 
in IRX between the two epochs. The right combination of grain sizes, gas density, and ion 
exposure time is therefore required to produce the observed dust temperature and decrease 
in IRX. 



4.4.2. The Grain Destruction Rate in the Hot Plasma 

Figure [9] depicts the evolution in the IRX{t2) / IRX{ti) ratio as a function of ti — to, 
where to is the time, since the explosion, when the SN blast wave first crashed into the 
dense material of the ER, taken here to be the independent variable. The times ti = 6190 d 
and ^2 = 7137 d correspond to the two epochs of near-simultaneous Spitzer and Chandra 
observations of the ER. The I RX {t2) / 1 RX (ti) ratio is derived by taking the the ratio of the 
fraction of the surviving dust mass, M^/M^, at the epochs ^2 — and ti — to (see Figure [6]). 
Implicit in the figure is a conversion from the dimensionless time variable t/Tmax to absolute 
time. 

The observed IRX{t2) / IRX{ti) ratio is ~ 0.53 ± 0.16 (see Table 1), and shown as a 
horizontal dashed line in the figure. The grain size distribution used in the calculations is 
characterized by a a~^'^ power law in grain radii, extending from a minimum grain size of 
10 A to a value of amax of 0.01, 0.04, and 0.1 /im. Results are presented for a SN blast wave 
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expanding into a one dimensional protrusion (a = 0, left column) and a uniform medium 
[a = 2, right column) with densities of 10'^ cm~^ (top row) and 10^ cm~^ (bottom row). The 
red line in the two bottom figures correspond to calculations performed for the grain size 
distribution of the ER with {amm, amax} = {0.023, 0.22} /im. 

The figure shows that the I RX (t2) / 1 RX (ti) ratio attains its lowest value when ti — to 
is small, that is, when the first encounter of the ER with the SN blast wave occured just 
before ti, the first epoch of Spitzer observations. Since ti is very close to Iq, very little 
grain destruction could have taken place during the ti — to epoch. The value of IRX{ti) 
is therefore close to its pre-shock value. Consequently, any subsequent destruction would 
lead to relatively great changes in IRX at t = t2- Conversely, IRX{t2) / IRX{ti) attains its 
largest value when to = (a physical impossibility because of the finite time required for the 
SN blast wave to reach the ER). Then, the relative change in IRX between the two epochs 
will be the smallest, and I RX {t2) / 1 RX (ti) 1. To illustrate the asymptotic behavior 
of the ratio of IRX values the curves were drawn for ti — to values beyond the maximum 
physical value of 6137 d. 

Figure [H] also shows the implicit dependence of the IRX{t2)/ IRX{ti) ratio on the grain 
destruction rate. When Tmax is small compared to the time scale of blast wave-ER interaction, 
the grain size distribution relaxes to its equilibrium form, and the curves of IRX{t2) / IRX{ti) 
for the different grain size distributions converge to the same functional form at sufficiently 
large values of ti — to. This is especially evident in the lower two figures, for which the 
gas density is higher, and grain destruction time scales are lower. At ti — to > 1000 d the 
curves of IRX{t2) / IRX{ti) for the different size distributions have all converged to the same 
functional form. 

The figure shows that the value of ti— to ranges from ^ 500—2000 d, for IRX[t2) / IRX(ti) = 
0.53 ± 0.16, giving values of to ~ 5700 — 4200 d for a = 0. The range of values for to is 
consistent with the uncertainties in the time of the blast wave encounter with the ring. This 
suggests that the grain size distribution derived from the X-ray constraint on the ionization 
timescale, the range of plasma densities, and the grain destruction efficiency used in the 
model are all consistent with the observed evolution of IRX, and estimated epoch of the 
first interaction of the blast wave with the ER. 

The problem can also be reversed to determine the grain destruction efficiency by ther- 
mal sputtering, by adopting the epoch of to = 4000 d and the nominal value of IRX{t2) / IRX{ti) = 
0.53 as accurate descriptions of the encounter time and the rate of decrease in IRX in the 
6190 to 7137 time interval. The the value of r^ax needs then to be adjusted to move the 
intercept from the calculated value of ti — to ~ 1000 d to the desired value of ~ 2000 d. This 
will require an increase in Tmax by a factor of 2, or a decrease in the destruction rate of the 
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silicate grains by a factor of ~ 2. 



4.5. Implications for Determining the Mass of the Circumstellar 
Environment of SN1987A Using Light Echoes 

It is interesting to compare the dust properties derived for the ER with those derived 
for the progenitor's circumstellar environment from studies of the evolution and intensity of 
light echoes created by the scattering of the optical light from the supernova by the dust 
grains. At any given time, all points with equal delay time lie on an ellipsoid of revolution 
with the SN at one focal point and the observer at the other. Unfortunately, the ellipsoid 
at the earliest epoch at w hich the echoes were observed was outside the ER [see Figure 11 



m 



Sugerman et al.l (j2005al )]. As a result, the light echoes probed only the circumstellar and 



interstellar media exterior to the ER. 

Assuming cylindrical and reflection symmetry, Sugerman et al. ( 2005bl lal) derived a 
model for the morphology of the scattering medium consisting of: (1) a peanut-shaped 
contact discontinuity (CD) between the red supergiant and main- sequence winds from the 
progenitor star; (2) a structure called Napoleon's Hat (NH) constituting the waist of this 
peanut; and (3) the two outer rings of the circumstellar shell (CS) that define th e hourglass 
that is pinched by the ER. T o model the scattered light ISugerman et al. J2005bl lal) used the 
Weingartner fc Draind (120011 ) model for interstellar LMC dust with grain radii ranging from 
an upper limit of 0.2-2.0 /im to a lower limit of 0.00035 fim. By varying the relative silicate- 
to-carbon dust mass ratio while maintaining an LMC dust-to-gas mass ratio that is 0.3 times 
the value of the local interstellar medium, they estimated a total nebular mass of 1.7 Mq. 
They also found that the gas density increases, the maximum grain size decreases, and the 
silicate-to-carbon dust mass ratio increases as the echo samples material that is closer to the 
SN. The ER, with its population of smaller pure silicate grains, is consistent with this trend. 
The higher value of the minimum grain size in the ER may be the result of its proximity to 
the SN which caused the evaporation of grains smaller than 0.02 /im by the initial UV flash. 



Finally , the dust abundance in the ER is consistent with that adopted by ISugerman et al. 



(j2005al ) for the nebula, supporting their derived value for the nebular mass. 



5. Summary 



The interaction of the SN 1987A blast wave with the complex structure of the ER has 
given rise to rapid evolutionary changes in the X-ray, optical and mid-IR morphology of the 
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emission. The Gemini South mid-IR images have estabhshed that the IR emission originates 
from dust in the ER that is swept up by the SN blast wave, and coUisionaly heated by a 
soft X-ray component which has a temperature of 3.5 x 10^ K, and an ionization timescale 
oi Uet > 10^ cm~^ d. The Spitzer infrared observations provide important complementary 
information on the evolution of the interaction of the SN blast wave with the ER and the 
properties of the dust in the hot X-ray emitting gas. The results of our analysis are as 
follows: 

1. Spitzer spectral observations on day 6190 after the explosion revealed that the dust 
consists of silicate dust grains radiating at an equilibrium temperature of ~ ISOi^g K. 
Subsequent observations on day 7137 revealed that the IR flux increased by a factor of 
~ 2, with the same dust composition and temperature remaining the same (Figure [7]). 

2. The narrow range of grain temperatures constrains the range of grain sizes and the 
combinations of plasma temperature and densities capable of heated the dust to the 
observed range of temperatures. Using the X-ray constraint on the ionization timescale 
we limit the grain size distribution in the preshocked gas to be between 0.023 and 
0.22 fim. The grain size distribution may have originally extended to smaller radii, but 
if so, these grains were evaporated by the initial UV flash from the SN. 

3. The observed value of IRX, the IR-to- X-ray flux ratio is consistent with that expected 
from a dusty plasma with LMC abundances of heavy elements. 

4. The value of IRX decreased by a factor of ~ 0.53 ± 0.16 between days 6190 and 
7137, suggesting that we are witnessing the effects of grain destruction on a dynamical 
timescale of the remnant. The magnitude of the decrease in IRX between the two 
epochs is strongly determined by the following factors: the grain size distribution in 
the preshocked gas of the ER; the density of the X-ray emitting gas and the grain 
destruction efficiency; and the total time the dust is exposed to the flux of sputtering 
ions. 

5. To follow the evolution of IRX, we constructed a model for the evolution of the grain 
size distribution in the shocked gas. In the model, pristine dust is continuously injected 
into the hot gas by the expanding SN blast wave, and destroyed by thermal and kinetic 
sputtering behind the shock. The evolution of the grain size distribution resulting 
from the combined effect of dust injection and destruction is presented in Figure [5] for 
different geometries of the medium into which the blast wave is expanding. 



6. 



The evolution in IRX represents the changes in the dust-to-gas mass ratio in the 
shocked gas resulting from grain destruction (see Figures [6] and [9]) . Given the grain 
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size distribution and plasma density, the decrease in IRX between two epochs can be 
used to determine the epoch at which the dust was first swept up by the SN blast 
wave. Conversely, knowledge of the epoch at which the SN blast wave first encounters 
the ER can be used to determine the grain destruction efficiency in the hot gas. 

7. A self-consistent picture that emerges from the application of the model to the com- 
bined X-ray and IR observations is that of a SN blast wave expanding into a dusty 
finger-like protrusion of the ER with a typical LMC dust-to-gas mass ratio. The dust 
in the preshocked gas consists of pure silicate dust with a normal LMC dust-to-gas 
mass ratio and a grain size distribution limited to radii between ~ 0.023 and 0.22 /xm, 
sufficiently large to stop the incident electrons. Smaller grain sizes may have formed in 
the mass outfiow from the progenitor star but were probably vaporized by the initial 
UV fiash from the SN. The SN blast wave crashed into the ER between days 4000 
and 6000 after the explosion giving rise to the observed soft X-ray emission. Typi- 
cal temperatures and densities of the soft X-ray emitting gas are ~ 3 x 10^ K and 
(0.3 — 1) X 10^ cm~^, consistent with those required to collisionally heat the dust to 
a temperature of ~ ISOit^g K. The plasma parameters and grain size distribution are 
consistent with the amount of grain destruction needed to account for the observed 
decrease in the IRX fiux ratio between days 6190 and 7137. At these gas densities, 
the onset of grain destruction occurred about 1200 - 2000 days before the first Spitzer 
observations, consistent with the rise in the soft X-ray fiux and the ionization time 
derived from X-ray models. 

Further Gemini, ESO VLT, and Spitzer observations of SNR 1987A are in progress 
which, with combined X-ray observations, will shed further light on the nature of the mor- 
phology and dust properties of the circumstellar medium around the SN. 
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with NASA. E.D. acknowledges partial support from HST grant GO-9114 for the Supernova 
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R.D.G. E.F.P., and C.E.W were supported by NASA through contract No. 1215746 issued 
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Table 1. Observed X-ray and Infrared Fluxes From SN 1987A^ 



day number^ 


X-ray fiux'^ 


fsoft 


IR flux 


IRX flux ratio^ 


6190 


(2.1 ±0.32) 10-12 


0.50 


(5.1 ±0.9)10-12 


4.9 ± 1.1 


7137 


(6.4 ±0.32) 10-12 


0.60 


(1.0 ±0.18) 10-11 


2.6 ±0.5 



iFluxes are in units of erg cm"2 g-i 
^Since the explosion. 

^X-ray flux in the 0.5-2.0 keV band, interpolated to the epochs of the Spitzer 
observations and corrected for an extinction column density of A'^^^ = 2.35 x 
10^1 cm"2. The error represents the uncertainty in the pile up correction factor 
(Park et al. 2007). 



4 



The fraction of the 0.5-2.0 keV flux that arises from the soft [kT ~ 0.3 keV) 



X-ray component (jPark et al.ll2005l ). 

^The ratio of the IR to soft X-ray flux from the SN. The value of IRX has 
decreased by a factor of 0.53 ± 0.16 from day 6190 to day 7137. 



